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The Enskog Process 


S. Albeverio*, B. Riidiger'^and P. Sundar ^ 


Abstract 

The existence of a weak solution to a McKean-Vlasov type stochastic differential system 
corresponding to the Enskog equation of the kinetic theory of gases is established under 
suitable hypotheses. The distribution of any solution to the system at each fixed time is 
shown to be unique. The existence of a probability density for the time-marginals of the 
velocity is verified in the case where the initial condition is Gaussian, and is shown to be 
the density of an invariant measure. 

Key WORDS: Boltzmann and Enskog equations; McKean-Vlasov stochastic differential 
system; invariant measure. 
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1 Introduction 


The Boltzmann equation describes the time evolution of the density function in a phase 
(position-velocity) space for a classical particle (molecule) under the influence of other particles 
in a diluted (or rarified) gas [lOj (evolving in vacuum for a given initial distribution). It forms 
the basis for the kinetic theory of gases, see, for e.g. ca¬ 
lf / is the density function, which depends on time t > 0, the space variable x G and 
the velocity variable u G of the point particle, then f{t,x,u)dxdu is by definition the 
probability for the particle to have position x in a volume element dx around x and velocity 
u within the volume element du around u. For a single type of particles all of mass m > 0, in 
the absence of external forces, the Boltzmann equation has the general form 

d f 

x,u)+u- X, u) = Qif, f){t, x, u), (1.1) 
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where Q is a certain quadratic form in /, called collision operator (or integral). 
Set A := X (0,7r] x [0, 27r). Then Q can be written in the general form 


Q{f,f){t,x,u)= f {f{t,x,v*)f{t,x,v*)-f{t,x,u)f{t,x,v)}B{u,dv,dO)d4>. (1.2) 

Ja 

It is assumed that any gas particle travels straight until an elastic collision occurs with another 
particle. Each u £ in (II.2h denotes the velocity of an incoming particle which may hit, 
at the fixed location x £ particles whose velocity is fixed as tt £ M^. Let u* £ and 
V* £ denote the resulting outgoing velocities corresponding to the incoming velocities u and 
V respectively. 0 £ (0, vr] denotes the azimuthal or colatitude angle of the deflected velocity, v* 
(see [E]). Having determined 0, the longitude angle (j) £ [0, 27r) measures in polar coordinates, 
the location of v*, and hence that of u*, as explained below. 


In the Boltzmann model as the collisions are assumed to be elastic, conservation of kinetic 
energy as well as momentum of the molecules holds, i.e. considering particles of mass m = 1, 
the following equalities hold: 


V* + v* = u + v 
*\2 


{u*y + {v*y = u^ + v 


(1.3) 


u* = u + (n, n — u)n 
u* = u — {n,u — v)n 


(1.4) 


where 


V* — V 

n = - -r 

|u* — v\ 

where (•, •) denotes the scalar product, and | • |, the Euclidean norm in 


(1.5) 


Remark 1.1. The Jacobian of the transformation (Q has determinant 1 and (u*)* = u 
since the collision dynamics are reversible. 


The outgoing velocity u* is then uniquely determined in terms of the colatitude angle 0 £ (0, vr] 
measured from the center, and longitude angle cj) £ [0, 27r) of the deflection vector n in the 
sphere with northpole u and southpole v centered at which are used in equation a and 
(Ha) (seee.g. [H], [21], |11|). 

B{u,dv,d9) is a a-additive positive measure defined on the Borel u-field B(M^) x B((0,7r]), 
depending (Lebesgue) measurably on u £ The form of B depends on the version of 

Boltzmann equation one has in mind. In Boltzmann’s original work uni, 


B{u, dv, dO) 


\{u — v) ■ n\dvd6 


u — v\dv cos 




de, 


( 1 . 6 ) 
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(1.7) 


where in (|1.6I) we used that f — f is the angle between u — v and n, so that 

|(M-'t;,n)| = |m-?;|cos(^ “ ^ W - v\ sm{^) , 

9 9 

and |(n — n, n)|(in = |u — n] sin(-) cos(-)(i0d(/> = il(u, dn, (1.8) 

is the differential cross section scattering the velocities v of incoming particles colliding with 
the particle with velocity u, written in polar coordinates. 

In the case where the molecules interact by a force which varies as the nth inverse power of 
the distance between their centers, one has 

n—5 

B{u,dv,d9) = |(n — t')|"-i/3(0)dt'd0 (1.9) 

where /3 is a Lebesgue measurable positive function of 6. In particular, for n = 5, one has the 
case of “Maxwellian molecules”, where 

B{u,dv,d9) = I3{9)dvd9. 

The function /3(d) decreases and behaves like for d | 0, see, for e.g. [H]. Note that 

in the latter case j3{9)d9 = +oo. We note that for Maxwellian particles the cross section 
B{u,dv,d9)d4> does not depend on the modulus \u — v\ of the velocity difference between the 
velocity u of the particle and the velocities v of incoming particles. 

In the present paper, we shall mainly assume that 

B{u, dv, d9) = a{\u — v\)dvQ{d9) (1-10) 

where Q is a cr-finite measure on ;S((0,7r]), and cr is a bounded, Lipschitz continuous, positive 
function on M"*". The assumptions on a, though an improvement on the existing results for the 
case (7 = 1, do restrict applicability to physically realizable molecules. When Q{d9) is taken 
to be integrable, one speaks of a cut-off function. 

Remark 1.2. Both the rigorous derivation of Boltzmann equation from a microscopic model, 
and the study of existence, uniqueness and properties of solutions of the Boltzmann equation 
still present many challenging and open problems. For the derivation problem, see, for e.g., 

m 

Morgenstern [32] “mollified” Q{f,f) by replacing it by 
QMif,f){t,x,u) 

= j {/(i, X, u*)f{t, y, V*) - f{t, X, u)f{t, y, v)}Km{x, y)B{u, dv, d9)dyd(j) 
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with some measurable Km and B such that Km{x., y)B{u, dv, dO) has a bounded density with 
respect to Lebesgue measure dv x dO, and obtaining a global existence theorem in x 

M^). Povzner [35] obtained existence and uniqueness in the space of Borel measures in x, 
with the term KM{x,y)B{u,dv,d9)dyd4’ replaced by Kp{x — y,u — v)dvdy (with a suitable 
reinterpretation of the relations between x, u* ,v* and y, u, v, and suitable moments assumptions 
on Kp). 

According to [13] (p. 399), this modification of Boltzmann’s equation by Povzner is close to 
physical realiti/^. Cercignani also notes that Povzner equation has a form similar to the Enskog 
equation for dense gases, which we shall discuss below. Modification in another direction 
consists in taking the space of velocities as discrete, and is discussed in [T3] (pp. 399-401). 

The majority of further mathematical results concerns the spatially homogeneous case, where 
the initial condition on / is assumed to be independent of the space variable x, so that at all 
times / itself does not depend on x. For such results see, e.g. m, [IS]. 

Let us now associate to mi its weak (in the functional analytic sense) version. The 

following proposition is instrumental in this direction. 

Using also R,emark ll.l[ Tanaka [32] proved the following result that is important for the weak 
formulation of the equation. 

Proposition 1.1. Let 'k(x,u) G (70(1^®), as a function of x £ u G M^. With B as in 
(fTTUD . we have 


/ ^(x, u)f{t, x, u*)f{t, x, v*)B{u, dv, d6)dxdud(l) 

Jm?x(0,tt]x[0,2-k) 

= / ^{x,u*)f{t,x,u)f{t,x,v)B{u,dv,d6)dxdudcl) 

4 ]R 9 x( 0 , 7 r]x[ 0 , 27 r) 


( 1 . 11 ) 


The above result is proven using equation dm) and Remark 11.11 |42j . From now on, we will 
assume that B is as in (|1.1U|) . 


Weak formulation of the Boltzmann equation 


Consider the Boltzmann equation with collision operator (II.2p . We multiply (dTD by 
a function (of {x,u) G M®) belonging to Cq(M®), and integrate with respect to x and u. 
Using integration by parts and Proposition we arrive at the weak form of the Boltzmann 

equation: 


f if{x,u)-^{t,x,u)dxdu — f f{t,x,u){u,Vx'f>{x,u))dxdu 
= / f{t,x,u)Lf'if{x,u)dxdu (1.12) 
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for all t G M_|_ with 


Lf'il;{x,u) = / —'ip{x,u)}f{t,x,v)B{u,dv,d9)d(l), 

X {0,7r] X [0,27r) 

where B is as in (|1.10l) . 

To proceed further, let us introduce an approximation to the weak form of the Boltzmann 
equation by introducing a smooth real-valued function (5 (which should not be confused with 
the one appearing in ()1.9I) 1 with compact support defined on 


f 'ip{x,u)-^{t,x,u)dxdu — f f{t,x,u){u,Vx'il^{x,u))dxdu 
Jrs at 

= [ f{t,x,u)L^Mx,u)dxdu 


(1.13) 


for all if} G (^(^(M®) and for all t G M+ with L^j:'ip(x,u) 

= [ {'ilj{x,u*)-'ip{x,u)}f{t,y,v)/3{\x-y\)dyB{u,dv,d9)d4>. 

J]R6x(0,7r]x[0,27r) 

Heuristically, when (3 ^ 6o, then any solution of (j 1.1311 tends to a solution of Boltzmann’s 
equation (|1.12p . so that fd can be seen as a regularization for (I1.12p . 

Equation (jl.l3p is thus the (functional analytic) weak form of an equation closely related to 
the Boltzmann equation, which can be written as 

x,u) + u- Vxf{t, X, u) = QeU, f){t, X, u), (1.14) 

with 


QE(/>/)(iia:,u) 


{f{t, y, X, V*) - f{t, y, u)f{t, X, v)}/d{\x - y\)dyB{u, dv, d9)d(j). 


In the case where jd is replaced by the characteristic function (or a smooth version of it like in 
m ) of a ball of radius e > 0, this is Enskog’s equation used for (moderately) “dense gases” 
taking into account interactions at distance e between molecules. For Enskog’s equation, see e.g. 
CD, US], m (pp- 6, 14 ), m, cs], la, a, m, m, a, m- For versions of the equation in a 
bounded region, see The relationship between the Enskog and the Boltzmann equation 

have been discussed in several publications. In particular, their asymptotic equivalence (with 
respect to the support of /d shrinking to {0}) has been discussed in [6]. In [37], a pointwise 
limit has been established. 
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If /it denotes the Borel probability measure on corresponding to a smooth density function 
f{t,x,u), i.e. 

Iitidx,du) = f{t,x,u)dxdu, 
then the equation (jl.ldp can be written as 

d 

V’) - (u, V^ijj{x, u))) = (/it, (1.15) 

where 

= / {ip{x,u*) —'ijj{x,u)}f3{\x — y\)nt{dy,dv)B{u,dv,d6)d(t>. 

Jr 6 X (0,7r] X [0,27r) 

In the above, we have used the sharp bracket (•, •) to denote integration with respect to /it 
while (•,•) denotes the inner product in If /it satisfies (ll.ldp . we say that /it is a weak 
solution of the Enskog equation. 

Define the space D := D(R+,M^) as the space of all right continuous functions with left limits 
defined on [0, oo) taking values in M^, and equipped with the topology induced by the Skorohod 
metric (see e.g. i). We denote the value of any cj S D at any time s by ujg or a;(s). Likewise, 
the time marginal of a Borel probability measure /i on B will be denoted by /is for all s G [0, oo). 
The measure /is will be a Borel probability measure on We will use similar notations for 
functions in B x B and for Borel measures on B x B. 

If the measure /if in (11.151) is the marginal at time t of a Borel probability measure /i on B x B, 
then we can write the Enskog equation (|1.15l) as follows: 

d 

— (/i, 'tp{xt,ut)) - (/i, (iif, Vcct'tp{xt,ut))) = (/i, Lf^'ilj{xt,ut)) (1.16) 

where xt,ut are t-coordinates in D x D and L^'il>{xt,ut) 

= / {'ijj{xt,u*) -'ip{xt,ut)}a{\ut - vt\)P{\xt - yt\)y{dy,dv)Q{d9)d(f), 

JUo 

where we used the form (ll.lOp of B and the notation Uq for B x B x (0,7r] x [0,27r). In the 
following, we formulate the connection between equation (|1.16l) and stochastic analysis. 

Remark 1.3. The idea of looking at solutions of certain deterministic nonlinear parabolic evo¬ 
lution equations in connection with probability measures describing the distributions of suitable 
associated Markov processes goes back to McKean fgP| /. For a spatial homogeneous version of 
our present context for the case <7 = 1, this idea has been adapted and ingeniously implemented 
by Tanaka and successively developed, for this case, e.g., in ® agf, m ,« 

l^ . In our work, we avoid the assumption of spacial homogeneity, and we allow a to depend 
on \u — v\. 
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Let us first derive heuristically the evolution of the stochastic process (Xs,Zs)sg]R, describing 
the evolution of position and velocity of a particle evolving according to the Enskog equation 
(jl.lGD . In the present context, the evolution of the velocity (Zs)seR of one particle is obtained 
by integrating (or in other words ’’summing”) the velocity displacements a{Zs,Vs-9,4>) with 
respect to a counting measure Nx,z{ds, dy, dv, d9, dcp) which depends over a time interval ds on 
the distribution y{dy, dx) in position and velocity of the gas particles, as well as the position and 
velocity (X^, Zg) of the particle itself because of the presence of particles being close enough to 
hit (guaranteed by the function /3), and the scattering measure for the velocity B{u, dv, d9)d(j), 
defined in (jl.lOh . The position then evolves according to Xt = Xq + Zgds. Let us introduce 
such a suitable jump-Markov process (dfs,Zs)sgR. Let fj,{dx,dv) be a probability measure on 
D X D. Let Nx,z{ds,dy,dv,d9,d(p) be a compensated random measure (crm) defined on a 
filtered probability space (n,^,^t,P) with compensated measure (or simply, compensator), 

dr := T{dy,dv,d9,d(p,ds) = a{\Zs — Vs\)l3{\Xs — ys\)fJ^idy,dv)Q{d9)d4>ds (1-17) 


on B X B X (0, tt] X [0, 27 r) x ]R_|_ (where we recall that Q{d9) is a cr - hnite measure on the Borel 
a - algebra ;8((0, vr]), dcj) is the Lebesgue measure on i3([0,27r)), and /3 is a C“(]R^) function 
with support near 0. Here, X, Z are elements of B, and Vs,ys are s-coordinates of v,y in B. 


Now, let {Xg, Zs)sgR^ be the process defined below, taking values in the Skorohod space B x B 
with the joint distribution of (X^, Zs)seR^ denoted by y{dx, dz): 


Zt —Zn + 


fl 
fl 


DxBx (0,7r] X [0,27r) 


a{Zg,Vg, 9, (p)Nx,zids, dy, dv, d9, d(p) 


+ 11 a{Zg,Vg,9,(j))a{\Zg - Vg\)P{\Xg - yg\) 

10 dlD)xDx(0,7r]x[0,27r) 


y{dy, dv)Q{d9)d(f)ds 


Xt —^0 + / Zgds 


/ 


The initial values Xq and Zq are specified. We have set 

a{u, V, 9, (p) := (n, u — v)n. 


(1.18) 

(1.19) 


( 1 . 20 ) 


where, as above, the deflection vector n is given in spherical coordinates, i.e. in terms of the 
colatitude angle 9 G (0,7r] and longitude angle (p G [0,27r). 

We have obtained such a process heuristically considering the physics governing the evolution 
of the particles and will prove in this article that this is the stochastic process whose law 
corresponds to the solution of the Enskog equation (|1.16l) . We remark however that the 
stochastic equation (II.1811 . (|1.191) is dehned in terms of a counting measure Nx,z with random 
compensator a{\Zs — Vg\)/3{\Xg — ys\)fj,{dy,dv)Q{d9)d(pds. The mathematical theory of point 
processes with random compensator has been analyzed extensively in e.g. [25], or 12ZI, but as 
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the theory of Stochastic Differential equations with Poisson random measure is more developed 
and better known, we prefer here to rewrite ()1.18|l . (ll.l9p in an equivalent stochastic equation 
written in terms of a stochastic integral w.r.t to a Poisson random measure associated to a 
Levy process, which is the following: 


t [ aiZs 

Jo Ju 


Zt = Zq + I I 
>0 Ju 

rt 


Vs, 0, (/>) 1 [ 0 , a{\Zs-Vs\)Ji\Xs-ys\)]iv)dN^ 


+ [ [ a{Zs,Vs,6,(l))(7{\Zs - Vs\)/3{\Xs - ys\)Kdv,dy)Q{de)d(l)ds 

Jo Juo 


( 1 . 21 ) 


Xt = Xo + 



( 1 . 22 ) 


where p is still the law of the process {Zt, Xt), t > 0, but now N^{dy, dv, dOdcf), dr, ds) is a com¬ 
pensated Poisson random measure (cPrm) with Poisson measure :=NV(^dy, dv, dOdcj), dr, ds) 
and compensator //(dy, dv)Q{d9)d4>drds on D x D x (0, tt] x [0, 27r) x [0,1] x M_|_ . 


That (ll.2ip . (ll.22p and (|1.18p . (ll.l9p are equivalent equations can be shown with at least two 
different methods: 


i) For each Borel-subset i? of B x D x (0, vr] x [0, 27r), the counting measure Nx,z{B x [0, t)) 
can be represented by 

[ [ [ Mo, ^{\Zs-vsm\Xs-ys\)]ir)dNV . 

Jo JbJ[o,i] 

This is a consequence of the following equation, which shows the relation between the 
random compensator T of the point measure Nx,z defined (I1.17p and the compensator 
fi{dy, dv)Q{d6)d(j)ds of the Poisson random measure N^: 


T{B X [0,t)) = [ f (t{\Zs -Vs\)I3{\Xs - ys\)iJ,{dy,dv)Q{d9)d4>ds 
Jo Jb 

^[o,ai\Zs-vs\)di\Xs-ys\)]i'^)d-idy,dv)Q{de)d())dsdr 



0 Jb j[o,i] 


ii) Given iIj G (^^(M®), for all t > 0, At > 0, the Ito formula ^p{Xt+At, Zt+At) — 'M{Xt,Zt) 
for (jl.2ip . (jl.22D and (jl.l8h . (ll.l9p can be proven to be exactly the same. This implies 
that p.21|) . (I1.22P and (I1.18p . (jl.l9l) solve the same martingale problem and are the same 
process in weak sense. In particular the law of the corresponding process {Zt,Xt), t > 0 
solves in both cases the Enskog equation, as proven in Proposition 12.21 for (ll.2ip . (ll.22l) . 
The Ito formula for (jl.2ip . (I1.22p is computed in the proof of Proposition 12.21 It can be 
obtained in a similar way for (I1.18p . (jl.l9l) by using Theorem 2.42, Ch. II in |27j . 












































We call the process {Zt,Xt), t > 0, given by (jl.21h . (ll.22p . (ll.2np lresD. its law fi = fit,t > 0) as 
the Markov process (resp. law) associated with the Enskog equation described by (|1.16p . Its 
existence is proven in Theorem l2.1l in Section 2 under suitable conditions which are satisfied by 
some physical models. In Proposition 12.21 we will prove that for any finite fixed time T > 0, its 
law fi = {/rt})0 <t<T solves the Enskog equation (II.IBp . Uniqueness of the Markov process 
{Zt,Xt), t > 0 solving (II.2ip . (II.22p . ()1.20p . is proven in Theorem 17. 1 1 in Section 3 for the time 
marginals. The existence of a density f{t,x,z) for the distribution /r = {/r4},0 < t < T is 
proven in Section 4, for the particular case where the velocity marginals are time invariant. 
f{t,x,z) solves then the Enskog equation ()1.13p 

It is worthwhile to mention that we have not made the assumption of space homogeneity. We 
allow a that appears as the differential cross section (see equation (jl.lUI) l to depend on \u — v\. 


2 Existence Results 


In this section we establish the existence of a solution of the system of stochastic equations 
(jl.2ip . (ll.22p . with denoted by ^ that takes values in the set H := (0, vr] x [0,27r). Also, 

Q{d6)d(p is written as Q{d^) for notational simplicity. Prom the physical model, we know that 
Q{d^) should be a u-finite measure, and hence taken as a-finite. 

Hypotheses A: 

Al. The measure Q is hnite outside any neighborhood of 0, and for all e > 0, Q satishes 

9Q{d9) < oo. 

A2. a : —)• (as entering (|1.101) 1 is a bounded, Lipschitz continuous function on R+. 

There are many useful consequences of Al. Recall that a{z,v,^) = {n ■ {z — u))n with f — | 
as the angle between the vectors {z — v) and n (see 0). Hence, condition Al implies that 
there exists a constant C such that the following estimates hold. 


J \a{z, v,() - a{z', v',()\‘^Q{d() < C{\z - z'\'^ + [u - v'\^) (2.1) 

\a{z,v,C) - ot{z',v',i)\Q{dCj < C{\z - z'\ + \v - u'|). (2.2) 

Prom (j2.ip . it follows, by setting z' and v' to be z, and using the fact that a{z, z,^) = 0 for all 
G E, z G R^ that 

Jja{z,v,C)\‘^Q{dC) <C\z-v\‘^, (2.3) 
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and hence 


\a{z,v,0fQ{d^) < C{\z\^ + |up). 

(2.4) 

In a similar way, from ()2.2p it follows 


\a{z,v,C)\Q{d^) < C\z - v\, 

(2.5) 

and hence 


/ \a{z,v,f,)\Q{dC) <C{\z\ + \v\), 

(2.6) 


Condition A2 on a is required for mathematical reasons. It is worthwhile to note that for 
Maxwellian molecules, a is the constant function identically equal to 1. Hence hypothesis A2 
leads to more mathematical generality but still falls short of physical reality. 

Before we proceed further, we recall the following: Since the function /3 that appears in (|1.13l) 
is held fixed and has been assumed to be bounded, we will set ||/3||oo = 1- We will also take 
11 1 loo = 1- Besides, we take the constant C that appears in the estimates (12.ip - (12.6p to be 
greater than 1 in order to avoid writing C V1 in many of the estimates in this paper. A generic 
constant will be denoted by K though it may vary from line to line. 


Let us fix a finite time T > 0, and denote the Skorohod space D([0, T];M^) by D. We consider 
it here equipped with the Skorohod topology. Given a probability measure /r on B x D, let 
denote its marginal at time t. We define 


a{z,v,i) := a{z,v,i)a{\z - u|) 


for all z,v G and ^ G H (the function a was defined in (ll.20p L The main result of this 
paper is stated below. 


Theorem 2.1. Suppose that a is in C^piR). Suppose hypothesis A hold. Let Aq and Zq be 
valued random variables with finite second moments. For any fixed T > 0, there exists a 
stochastic basis (12,^, (=^t)iG[o,T]! P)? adapted process {Xt, Zt)t^yQ^T] values on B x B, 
and a compensated random measure (crm) N^, with p being the law of the stochastic process 
{X,Z), satisfying a.s. the following stochastic equation for t G [0,T].- 


Zt — Zq + 


fl 


0 JDxID)xHx[0,l] 


a(Zs, Vs, a{\Zs-vs\)h(\Xs-ys\)] {r)dN^ 


+ [ [ a{Zs,Vs,^)fii\Xs-ys\)dpQ{d^)ds 

Jo JDxDxS 


(2.7) 


and 


Xt = Xo + 



( 2 . 8 ) 


where dN^ :=N^{dy,dv,df,,dr,ds). For any t G [0,T], Xt and Zt have finite second moments. 
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Proposition 2.2. Let fi denote the law of the process {Xs,Zs : 0 < s < T}, solving ()2.7p . 
(j2.8l) . Then pL solves the Enskog equation M.lb]) for any € Cq(]R®). 


Proof. As for any t £ [0, T] Xt and Zt have finite second moments, and due to the conditions 
(Al) and (12.4p we can apply the Ito formula to (A*, In fact let t, At > 0, then 


V’(A^i+At) Zt+At) 

/ t+At 

{Z,,V^i^{Xs,Zs))ds 

rAt 


+ [ [ {ip{Xs,Zs + a{Zs,Vs,f)l[0^a{\z,-v,\)/3{\Xs-ys\)]ir))-'ll^{Xs,Zs) 

Jt Jc/ox[0,l] 


- V^fj{Xs, Zs)a{Zs,Vs,0, (/>)1[0, a(\Zs-vs\)gi\Xs-ys\}]i'^)}Lidy, dv)Q{d9)d(j)dsdr + 

[ [ a{Zs,Vs,9,f)),Vz'il>{Xs,Zs))}cr{\Zs -Us|)/3(|As - ys\)n{dy,dv)Q{d9)d4>ds 

Jt JUo 

ft+At 

t.S 

It 

+ [ [ {'ijj{Xs,Zs + a{Zs,Vs,f,)) -'iltiXs,Zs) 

Jt JUo 

- XzftiXs, Zs)a{Zs,Vs,9,f))}a{\Zs - - ys\)fi{dy, dv)Q{d9)d4>dsdr + 

rAt 


/ t+At 

{Zs,V^ft{Xs,Zs))ds 


[ a{Zs,Vs,9,(l)),X:,'ip{Xs, Zs))}a{\Zs -Us|)/3(|As - ys\)pi{dy,dv)Q{d9)d(j)ds 
JUo 

+ (2.9) 

where is a martingale, for each T G M. 

Taking the expectation E with respect to the measure fi{dx,dv) in ()2.9p . we get 


E[tlj{Xt+At,Zt+At)-ip{Xt,Zt)] 
rt+At 


+ 


+ 


/ r+ZAt 

iZsyMXs,Zs))ds] 

E[f f a{\Zs - ud)/3(|As - ys\){a{Zs,Vs,9,4>),V:,'if{Xs, Zs))pL{dy,dv)Q{d9)d(l)ds] 
Jo JUo 

E[f [ {f;{Xs,Zs + a{Zs,Vs,9,(P))-tP{Xs,Zs) 

Jo JUo 


- {Xz'if{Xs,Zs),a{Zs,Vs,9,(f))}a{\Zs - Us|)/3(|As - ys\)p,{dy, dv)Q{d9)d(j)ds] 


Dividing by At and letting At —)• 0 we obtain (|1.16l) by noting that /r is also the law of 
{Z,X). □ 
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In Section 7 in Theorem O we will prove uniqueness of the law ^ of the process {Xg, : 0 < 
s < T}, solving the McKean -Vlasov equation (12.71) . (12.8p in the following sense: we prove that 
for any fixed t in the interval [0, T], the t-marginal distribution of weak solutions of (12.71) . (|2.8I) 
is unique within the class of Borel probability measures on that are absolutely continuous 
with respect to the Lebesgue measure on M®. 

The existence of a probability density for the time-marginals of the velocity is verified in the 
case where the initial condition is Gaussian, and is shown to be the density of an invariant 
measure in Section 8, Theorem 18.II 


3 Existence and uniqueness of a stochastic equation 


Consider a given filtered probability space (fl, ^, (=^t)tG[o,T]) P) satisfying the usual conditions. 
Let St '■= denote the linear space of all adapted cadlag processes with values 

on equipped with norm 

ll^llsi sup |V,|]. (3.1) 

sG[0,T] 

is a Banach space. This can be shown similar to the proof of Lemma 4.2.1, page 93 in 
m- We consider for d = 6. 


Let B be equipped with the Skorohod topology, where B denotes B([0, r];M^). Given a prob¬ 
ability measure A on B x B, let At denote its marginal at time t. Let us assume 

/ [ i\vt\ + \yt\)Hdv,dy)dt < oo Vr>0. (3.2) 

Jo JdxO 

Consider a Poisson random measure on (11, ^tj P) with intensity measure Xt{dy, dv)Q{d^)drdt 
on the Borel subsets of B x B x (0, tt] x [0, 27r) x [0,1] x [0, T]. We denote by the corre¬ 
sponding compensated Poisson random measure. From the condition (|3.2p and (j2.6p . and the 
hypothesis ||cr||cxD= ||/3 ||cxd = 1, it follows that 



a{z,vt,C)W{\z - vt\)/3{\x 


yt\)X{dv,dy)dtQ{d^) < oo, 


(3.3) 


Vz S G M^,Vr > 0. 


Let us use the following notation: 


where B^ := B x B. 


[/o = B^X“ 

1/ = B^ X H X [0,1] 
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Theorem 3.1. Let {Zq,Xq) be a random vector with values on with 


E[|Zo|]<oo, E[|Xo|]<oo, 


(3.4) 


and assume Then for all T > 0 there exists a unique strong solution of the stochastic 

equation 


= Zo 

-I i 

+ [ a{Z^,Vs,0l3i\^s - ys\)Hdydv)Q{dC)ds 

JUo 

= Xo+ [ Z^ds, 

Jo 

on 5^, where dN^ denotes N^{dy,dv,df,,dr,ds). 


(3.5) 

(3.6) 


We first introduce some notation and preliminary results. 

Let T > 0 and (Z, X)ig[o^ 7 ’] be an adapted process with values in B^. 

Lemma 3.2. Assume (Z, X)jg[o,T] £ ^t- Then the stochastic integrals (L(Z))jg[Q and 
(.Hz))te[o,T], with 

I{Z)t:= f [ a{Zs_,Vs,f)l[o,a{\z,_-vs\)y{\x,-ys\)]{r)dN^, (3.7) 

Jo Ju 

and ^ 

i{Z)t:= f f \a{Zs_,Vs,f,)\l[o^a{\Zs_-vs\)y{\Xs-ys\)]i'^)dJ^^, (3.8) 

Jo Ju 

are well defined, and there exist constants K > 0 and Mt > 0 satisfying 

E[J(Z))t] < iL / E[sup |Z,|](it + MT. (3.9) 

Jo sS[0,t] 

Proof. We need to prove only inequality (13.9p . It then follows that the stochastic integrals 
(/(Z)))jg[o,T] and {I{Z)))t£[o,T] are well defined (see Section 3.5, in particular Lemma 3.5.20 
of [3^, or Theorem 4.12 [39]). 
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Using ()2.6h it follows 


E[/(Z)r] = / / E[\a{Zs,Vs,C)W{\Zs_ -Vs\)/3{\X^ -ys\)]X{dydv)Q{dC)ds 

Jo JUo 

<c( f E[|Zs|](is+ f [ \vs\\{dv,dy)ds\ 

\Jo Jo JoxB / 

<c( f E[ sup \Zs\]dt + [ f \vs\X{dv,dy)ds 
\Jo s€[0,£l Jo ^DxB 


XS < OO 


□ 


Let (5Z, 5X)jg[o^'r] denote the process defined through 

SZt := Zo + I{Z)t and SXt := Xq + [ Z^ds. 

Jo 

If the random vector (Zo,Xo) satisfies (13.4h and {Z,X) S then {SZ,SX) £ S^. 

Indeed 

sup \SZt\ < i{Z)T+ \Zq\ a.s. (3.10) 

te[o,T] 

and hence 

E[ sup \SZtW < E[/(Z)t] + E[|Zo|] (3.11) 

tG[ 0 ,T] 

The statement follows from the estimate ()3.9p . 


Lemma 3.3. For any T > 0 fixed, there exists a constant K > 0, such that for a// n £ N and 
all (Z,X) £ satisfying supgg[o,T] \^s\ < n, sup^gjo,T] \^'s\ ^ following inequality holds 

rt r 


0 JU 


\a{zs,Vs,i)l[Q^ a{\Zs-Vs\)y([\Xs-yA\i^) ~ ff(|Z'-i;.|)/ 3 (|X'-y«|)] (^) I 

X{dv, dy)Q{dff)drds 


< 


f Ln{\Zs-Z',\ + \Xs-X',\)ds P-a.s., 

Jo 


with Ln = Kn. 


Proof. 


n W{Zs,Vs,fi)l[Q^ a{\Zs-Vs\)y{[\X,-ysW{r) “ “(^s> , 01 [0, u{\Z',-Vs\)^(\X',-ys\)]{r)\ 

r 

X{dv,dy)Q{df)drds < I + II P — a.s., with 
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I : = 


0 JU 


a{Zs,Vs ,0 - OL{Z'^,Vs,i)\l[Q, a{\Zs-vs\)mXs-ysW{'^)Kdv,dy)Q{di)drds 

fc\Zs-Z',\ds 

Jo 


< 

where the last inequality follows from (12.21) . and 

rt 


II : = 


■■= f f \a{Z'^,Vs, 0 \ 

Jo Ju 


X 


{1[0, a{\z,-v,m\Xs-ys\)]ir) - Mo,a{\zi-v,\)fi{\x',-ys\)]{r)\}Hdv,dy)Q{d^)drds 

< [ [ \a{Z'^,Vs,0\\^^^i(^i\Zs-Vs\)l3{\Xs-ys\),cr{\Z'^-Vs\)Pi\X'^-ys\)) 

Jo JUo 

- va:ui{a{\Zs - ns|)/3(|Xs - ys\),a{\Z'^ - Vs\)/3{\X'^ - ys\)\)X{dv,dy)Q{dC)ds = 

[ [ \ct{Zs,Vs,0\ \cri\Zs - Vs\)(3{\Xs - y^l) - cr(|Z' - Vs\)(3(\X'^ - ys\)\X{dv,dy)Q{dC)ds 

Jo JUo 

Using that a and P are Lipschitz continuous functions bounded by 1, as well as we get 

that there exists a constant K > 0, such that 


II < 


[ [ \a{Z's,Vs,^)\ 

Jo JUo 


X (||Z' - Usl - \Zs - fsll + ||X' - ys\ - \Xs - ys\\)X{dv,dy)Q{d^)ds 


< K 


f [ \a{Zi,Vs,m\Zi - Zs\ + \Xi - Xs\)X{dv,dy)Q{dOds 
Jo JUr^ 


10 JUo 

<K [ [ {\Zi\ + \vs\)\i\Z',-Zs\ + \Xi-Xs\)X{dv,dy)ds 
Jo Jo^ 


□ 


In the next Lemma we will use the local Lipschitz condition stated in Lemma 13.31 to prove a 
local contraction property of 5 on 5^. 

Lemma 3.4. For each n G N there exists constant > 0, such that 

E[ sup \SZs - SZ'^\] <Ln [ E[ sup {\Z,, -Z'^,\ + - X'^,\}]ds 

se[o,t] Jo s'e[o,s] 

V(Z,X)^g[o_r] , (^',^')sG[o,T] e S^, with sup \Zs\<n, sup \Z'J<n. 

se[o,T] sg[o,t] 
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Proof. 


E[ sup \SZ, - SZ',\] < 

sG[0,t] 


E[sup / / \a{Z,,_,Vs',OMo,-{\z^> 

sG[0,t] JO JU 

- a{Z'^,_ , Vs' , 0 l[o, -j/,,I)] (^) I 

<E[/ f \a{Z,'_,Vs',OMo,<^{\z,'_-vM\Xs'-ys'\)](^^ 
0 ^ u 

- a{Z'^,_ , Vs' , 0 1[0, ,7(|Z',_-J/,, I)] (^) I 

= E[/ f \a{Zs',Vs',i)llo,a{\Z^,-vAmX^,-yM^'^) 
Jo Ju 

- a{Z'^,,Vs',OMo, <t(|z',-v|)/3(|x',-j/,,|)] 


<Ln [ n\Z's' -Zs'\ + |X', - Xs'Wds' 

Jo 

<Ln [ E[ sup {|Z', - Zs'\ + |X', - Xs'\}]ds 
Jo s'e[o,s] 

where we have used Lemma 13.31 □ 


In the proof of the next theorem we will use the local contraction property in Lemma 13.41 to 
prove existence and uniqueness of a modification of the stochastic equation defined through 
(13.5D . (13.61) . The modified stochastic equation satisfies global growth and Lipschitz conditions. 

Let j G N, Bj := {z G : \z\ < j} and 


aj{z,v,C) 


1 + d{z, Bj) 


(3.12) 


where d{z,Bj) denotes the distance of G from Bj. 

Theorem 3.5. Let the random vector (Zq,Xq) satisfy (|3.4I) . For all T > 0 there exists a 
unique solution on of the stochastic equation 



X ho,a[\zff-v4mxf’^-y4)](^)^^(^y^dv,dC,dr,ds) (3.13) 

x^’^ =Xo+ [ Z^dds. (3.14) 

Jo 
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It then follows directly the statement of the following corollary: 

Corollary 3.6. Let the random vector {Zq,Xq) satisfy (13.41) . For all T > 0 there exists a 
unique solution on of the stochastic equation 




= Zo 


+ 

+ 


Jo <7i\zf^-Vs\)/3{\xfX 

[ a,{Z^’^,Vs,0^i\Z^’^ - VsimiXt^ - 

JUo 

Xo+ [ Z^’^ds. 

Jo 


ys\)dX{dydv)Q{df,)ds 


(3.15) 

(3.16) 


Proof. If the stochastic integrals in the stochastic equation (13.131) , (|3.14p are well-defined then 
the stochastic equation (I3.13D . (13.141) is equivalent to the stochastic equation (I3.15p . (13.161) . 
(See e.g. Chapter 5 [30]). □ 


Proof of Theorem 13.51 


Proof. We start by remarking that 


1 -I- d{z, Bj) 


and there exists a constant Kj >0, such that 


< min(j, l^l) 


l + d{z,Bj) l + d{z',Bj) ^ 


(3.17) 


Assume (Z, A) 4 g[o^ 7 ’] S S)p. As \aj{z,v,f)\ < \a{z,v,f)\ it follows from Lemma [32] that the 
stochastic integrals (7(Z))ig[o^7’] (7(Z))igjo,r]5 with 

fj{Z)t:= [ [ \aj{Zs_,Vs,O\'^[0,a{\Zs_-vs\)mXs-ys\)]i^)dN^^ 

Jo Ju 


and 


Ij{Z))t:= / Oij{Zs_,Vs,C)l[o,ai\Zs_-vs\)h(.\Xs-ys\)]i'<^)dN 

Jo Ju 


are well defined. Moreover, using {Ij{Z))t < {I{Z))t Vt £ [0,T] and (|3.9p . it follows that 

rT 


E[/j(Z))t] <k[ E[ sup \Zs\]dt + Mt , 

Jo s€[0,t] 


(3.18) 
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with K > 0, Mt > 0. Let {SjZ, SjX)^^[Q xj denote the process defined through 

SjZt := Zo + Ij{Z)t and SjXt := Xq + [ Z^ds . 

Jo 

If the random vector {Zq,Xq) satisfies (13.41) . then 

E[ sup \SjZt\] < E[Ij{Z)T]+n\Zo\] 

t€[0,T] 

and, due to the growth condition (I3.18p . {SjZ, SjX)i^[Q rp] G S^. 

From ()3.17l) and Lemma 13.31 it follows that there is a constant Lj > 0 such that 

n \(y.j{Zi Us, ^)1[Q^ a(\z—V3\)fi(\\x—yt\] (^) “ ^^[O, a{\z'—Vs\)l3{\x—y3\)] (^)l 

r 

X{dv, dy)Q{d^)drds 

< Lj{\z — z'\ + \x — x'\) for all z,z',x,x'G M^, 

Similar to Lemma 13.41 it can then be proven that the following inequality holds: 

E[ sup \SjZs - SjZ'.W < Lj [ E[ sup {\Zs> - X',| + |Xs^ - X',,\}]ds 

sgfO.4 Jo s'gfO.sl 


(3.19) 


It follows that there exists n G N such that {S^Z, X)i^[Q rp] is a contraction from to S^. 
It then follows, that the mapping Sj has a unique fixed point on S^. □ 

Proof of Theorem 13.11 


Proof. It is sufficient to prove that for all T > 0 there exists a unique process {Z, X)ig[o,r]G (Sy 
satisfying a.s. the following stochastic equation 

= Zo+ 

IL a{Z^, Ol[ 0 , a{\z^-vsm\x^-ys\)]i'^)J^^idy, dv, dC, dr, ds) (3.20) 

X^ = Xo+ [ Z^ds, (3.21) 

Jo 

(13.201) . (13.211) is then equivalent to (13.5D . (13.61) . 

We will follow the strategy of the proof of Theorem 4.3.1 in |3D]. Let {Z^d S'!, 
be the unique solution of (13.131) , (j3.14p . Let 

Tj := inf{t G [0,T] : \zl’^\ > j] 
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By uniqueness of the solution of (j3.1,Sp . (|3.14p it follows that 


a.s., for t G [0,rj] , 

giving P(rj < Tj+i) = 1 Vj G N . 

We will prove 

P(UjeN{T,- = T]) = 1. (3.22) 

It then follows that the a.s. limit process {Z^,X^)f(z\nT^= hmi^cxDf-Z’'''’-^,is the 

solution of 630, eiD, and hence lEO, 

It follows from (j3.I9h and (I3.18P that 


]E[ sup \Z^'^\] <K [ E[ sup \Z^'^\]dt + Mt + E[|Zo|] 
tG[ 0 ,T] Jo se[ 0 ,f] 

so that by Gronwall’s Lemma 


E[ sup \Z^'^\] < expKT{MT + E[|Zo|]) 
te[o,T] 


(3.23) 


(3.24) 


It follows 


P{Tj < T) = P( sup \Z^’^ > j) 

tG[0,T] 

< -E[ sup \Z^’^\] < -expiLr(Mr + E[|Zo|]) 

J iG[0,T] J 

SO that 

P(GjGN{'L' < ^}) = 1™ < ^)} = 0- (3.25) 

j->-CXD 

□ 


4 Tightness 

In this section, we formulate an approximating sequence {Z^'^\ for the McKean-Vlasov 

limit, and prove the tightness of this sequence by Kurtz’s criterion on the Skorohod space B x B 
with the Skorohod topology. 

Define the processes 


Zf) = Zo 

XP =Xo + Zot 
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for all t G [0,T]. Let := Consider a Poisson random measure iV^(o) on 

U X [0,T] whose compensator measure is given by d^^^'>Q{d^)drds. Let -/V^{o) denote the 
corresponding compensated Poisson random measure (cPrm). By the square integrability of 
Zq and Xq, one has 


E( sup 

0<t<T L 


1 ^. 


( 0)|2 


+ 1 ^, 


( 0 )| 2 ' 


< oo. 


This implies in particular that VT > 0 {Z^^\X^^^) G S'^ C S^, where denotes here the 
Banach space of all adapted cadlag processes {Xt)t£[o,T] with values on equipped with norm 


XII 52 := (E[ sup |X,|2])V2 

sG[0,T] 


(4.1) 


(see e.g. page 93 in [3D]). In particular it implies that the measure is square integrable, 
and as a consequence satisfies the assumption ()3.2p . For all n > 0, define 

+ f [ -VsDdfx^^^QidOds (4.2) 

Jo JUo 

and 

X{n+i) =Xo+ [ (4.3) 

Jo 

Here, is the law of {Z^^\ and N^(n) is the cPrm with compensator measure given by 

d^jJ'^'^Q{d^)drds. Taking n = 0, it follows from Theorem 13.II that there exists a unique solution 
of {Z^^\X^^'>) in 5^ solving (14.2p . (14.31) . Moreover, 


E(|Z, 
< 3E 


W| 2 ^ 


+ 


< 3E 


(|Zop)+/ [ |a(zW,u„0lMl^i'^-^s|)/3(|XP-y,|)V)Q(d0d^ 

Jo JUo 

f [ a{Z(^\vs,OU{\Xi^^ -ys\)dfx<^°^Q{dOda 
Jo JUo 

(|Zo|2)+ f [ |a(zW,n„0lV°^Q(^^0rf^ 

Jo JUo 

t f [ \ [ a(zW,u„ 0 Q(^^ 0 PV°^ [ k(|zW-n,|)/3(|xW-2/,|)|2V0)d. 
Jo Jo^ Je Jo^ 


is T 
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by Cauchy-Schwarz inequality; continuing, by a use of (I2.1jl . (12.2p . and ()2.4h . one obtains 


< 3E 

< 3 


{\Zq\^) + C f[\zP\^ + \vs\^W^^Us + Ct + 

Jo Jo 


E|Zop + C(l + r) / E|zf)| 2 ds + Ct(l + r)E|Zo|^ 
io 


Hence by the Gronwall inequality, we obtain 

n\z[^'^?)<KiE{\Z^\‘^){l + K2t) 


(4.4) 


with Ki = 3e3C'r(i+G and K 2 = C{1 + T). It then follows for n = 1 that the law satisfies 
the assumption ()3.2p . so that there exists a unique strong solution {Z^^\ solving (14.2p . 

(14.3p . Along similar lines, one obtains 


E(|Z, 


< 3E 


( 2 )| 2 ^ 


l^or + 


C [\\ZP\^ + i\zi^^f))ds + CT [\\zP\^ + i\zi^^\^))ds 
Jo Jo 


so that by the Gronwall inequality. 


E\zi^^\‘^ < KiE{\Zo\^){l + KiK2t + 


(4.5) 


Further iterations result that (|4.2p . (|4.3I) has a unique strong solution (Z^”), A^”)) and the 
bound 


i=0 


so that for all n G N, we have 




(4.6) 


|_^(^)|2 _|_ |;^ W |2 


By the definition of we obtain an upper bound uniformly in n for E 

This is uniform boundedness of the sequence at each fixed t G [0, T]. Using the Burkholder- 
Davis-Gundy inequality, and proceeding exactly as above, one obtains an upper bound K 
uniformly in n for E supo<f< 7 ’(|Z^^”^ p + \x['^\'^) . It follows in particular that {Z^^\ 

G C S^. As S^, p = 1,2 are not separable Banach spaces, tightness has however to be 
proven on the Skorohod space B x D with the Skorohod topology: 
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in order to verify the second requirement in Kurtz’s criterion, we consider for any fixed 5 > 0, 
E 


I yi”-) _ yi”-) |2 I 77 


^t+s 

rt+S 


< 2E 


< 2E 


+ 


rt+o n 

rt+S r 

+ 1 / / d(zW,n„0/3(|xW-y,|)Q(dOV""'^o!sP}|-^i 

Jt JUo 

rt+S r 

{| / / la(zW,^;„0|V(|zW-r;,|)/3(|xW-y,|)V-i)Q(d0d5 

Jt JUo 

6 r' [ I [ a(z(”\r;„Og(dOl'V"'^x 

Jt Jd2 Jh 

[ a2(|zW-r;,|)/?2(|xW-y,|)V"-i)ds}|J-, . 

JD 2 

We will call the above expression on the right side as 2E(A^”^ | Tt)- Then, 

rt-\-5 

E(4''^ I J-f) < 2C(1 + 6 ) E(|zW| 2 + |zi”-i)|2)ds 

< K5 


for a suitable constant K > 0 which is independent of n. Hence, 

lim limsupE(A^"'^) = 0. (4-7) 

5^0 n^oo 

From (14.61) and (|4.7p . we conclude that is tight in B. By the definition of X^'^\ it follows 

that {Z^^\ is tight in B^. 


5 Distance Between Successive Approximations 

In this section, we give a result on the closeness of the measures and as n —>■ 00 . 

Proposition 5.1. Suppose that a is in C'^(M). Then for all h G C^piR®) and for any fixed 
t G [0, T], we have 

- h{zi''\xi'''^)] -+ 0 (5.1) 

as n ^ 00 . 

Proof. By the Ito formula, we write 'K[h{z[^~^^\x['^~^^'^) — h{z['^\x['^'^)] as Ai + A 2 where 

Ai=E[[ - V,,/i(zW,xW) • zWjds]. 

Jo 


22 








and 


A2=n / {h{Z(^+^^+a{Z(^+^\vs,C),X(^+^^)-h{Z(^+^\xi^+^^)} 
Jo JUo 

(’^■+ 1 ) _,, nfln t("'+i) - 


X a(|Zi-+i) - ^;,|)/3(|Xi-+i) - ys\)d^i^^^Q{dOds] 


Jo JUo 

X - 7;.1)/3(|xW - y,|)d/.(-i)Q(d0d5] 

By tightness, given any e > 0, there exists an i? > 0 such that 

P{ sup I inaxIlZp-*! + : j = n — 1,n,n + 1} > ii} < e. 

0<t<T 

This is a statement about the measures and Let Br denote the ii-ball in 

M®. First, we will deal with Ai. Clearly, by Cauchy-Schwarz inequality, 

I r f - f [ {V^/l(4”^4”^) 

Jo JJo Jb<S 


— Il^x^lloo 

<KtV~e 


JBg Jo 


1/2 


for a suitable constant Kt > 0. Restricted to Br, note that the function g{z, x) = Vxh{z, x) ■ z 
can be uniformly approximated by functions in C“(M®), so that 

Ai<Ky/e+ f sup I f — dfj,J^^}\ds. (5.2) 

Jo 0 gC^(k®) JK® 

In order to bound A 2 , we will split A 2 into two parts so that A 2 < Ii +12 where 

h=E[f [ {/i(zi-+i)+a(zi-+i),u„e),Xi-+i))-/i(zi-+i),xi”+i))} 

Jo JUo 

-E[f [ 

Jo JUo 

X a(|zW - U,|)/3(|XW - y,\)dfj,^^^QidOds], 
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and 


l2 = nf [ {h{Z^^^+aizi^\v,,^),xi^^)-h{zi^\xi^^)} 

Jo JUo 

X a(|zW - z;.|)/3(|xW - y^Ddfi^^^Q^ds] 

-nf [ {hiZ(-^+aiZ(-\v,,C),Xi^^)-h{Z^^\xi^^)} 

Jo JUo 

X a(|zW - z;.|)/3(|xW - y,\)dfi(^-^^Q{dOds]. 

To bound Ii, we will use the notation z,x instead of Zs,Xs in defining the function 


ip(z,x)= [ {h{z+ a{z,v,0,x) - h{z,x)}a{\z-v\)l3{\x - y\)dn^J^^Q{dC) 
Jr^xE 

where s is fixed. Given e > 0, there exists an i? > 0 such that for all n, 




e. 


A similar statement holds for in the place of Z^'^\ Therefore, we obtain 

[ |V^(z,x)|(d/rW+^i"+i))< Ae 

where AT is a suitable constant, and Bn is the i?-ball in M®. 


(5.3) 


Hence, we can focus our attention on Ii when the processes are restricted to values in Bn- 
Next, let Es denote the subset (0, (5] x [0, tt) of H. Given e > 0, there exists a 6 > 0 such that 
for all n and s, 

\[ [ {Hz + aiz,v,C),x) - h{z,x)}a{\z - v\)/3{\x - y\)fi^J^\dv,dy)Q{dC)\ 

JJ Bji X S5 


< e. 


(5.4) 


With this estimate in hand, we observe that the function of {z, x) given by 

[ {Hz + a(z, V, - Hz, x)}a(\z - v\)f3{\x - yDdHHQHH 

J Bji X (^,27r] X [0,7r) 


is a function in 


This observation along with (|5.3p and (|5.4I) yields 

|.^i|<A'e+ f sup I f (j){z,x){dy,H~^^^ — dyiHH\ds. 
Jo <AeC?°(K6) 4 k6 


(5.5) 


24 


For bounding I 2 , we repeat a procedure similar to the one used for Ji for the function 
giv,y) = [ {h{z + a{z,v,^),x) - h{z,x)}a{\z - u|)/3(|x - y\)Q{d^)^i^J^'>{dz,dx) 

JR®xS 

where s is fixed and the notation v,y is used instead of Vs,ys- We obtain the estimate 

\h\<Ke+ [ sup \[ (5.6) 

do 0gC^(IR®) jr® 

for a suitable constant K > 0. Combining (|5.2I) . (|5.5p . (j5.6p . we conclude that for suitable 
constants K, and Ci, 

\f h{z,x){dfi['"^^'^ - dn["''^}\ 

JR® 

<K€ + Ci f sup I f 4'{z,x){dfj,^J^~^^'^ — dfi^J^^}\ds 
Jo (AgC^(R®) JR® 

+ / sup I / 

JO 0gC^(R®) JR® 


We can take supremum on the left side over h € and call the resulting expression as 

Jn+i{t). Then we have 

Jn+i{t) ^ Ke + Cl f Jn+i{s)ds + f Jn{s)ds. 

Jo Jo 

By the Gronwall inequality, 

Jn+i{t) < iiTee^i* + [ Jn{s)e-^^^ds 


10 


I g-Cis 


< + C'le'-i" [ 

Jo 

< ... 


/ Jn-iir)e-^^^dr 

Jo 




rt rsi rsn-i 

+ / ••• / Ji{r)drds 

Jo Jo Jo 


1-1 • • - dsi 


/o Jo Jo 
which tends to zero as n —)• 00 and e —)• 0. 


□ 


6 Identification of the Limit 


In this section, we will conclude the proof of Theorem 12.11 by using tightness. Proposition 15.11 
and convergence of martingale problems. From the tightness of we have the existence 
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of a weakly convergent subsequence Let its weak limit be denoted by fi. Consider the 

associated subsequence which is also tight so that there exists a further subsequence 

which converges weakly. We will call its limit as u. Clearly, {//"'Ll being a subse¬ 
quence of converges weakly to fj,. Our aim in this section is to identify /r as a weak 

solution of the Enskog equation. 

Let us denote a generic element of B x B as cji x ^2. In the canonical setup on the path space 
B X B, we recall that is the solution of the following martingale problem: 

For any function (j) G x M^), and any t G [ 0 ,T], 


(i) ^ = C{Zo,Xo) (specified). 


(n) (p{uJi{t),LV2{t)) - (I){uii{0),lv 2{0)) - Vx(l){uii{s),u}2{s)) ■ UJl{s)d 6 

Jo 

n {())(a;i(s) +a{uJi{s),v{s),C),uj2is)) - (p{uji{s), u}2is))} 

^0 

(t(|wi(s) - Us|)/3(|w2(s) - ys\)Q{d()n'^''^i\dv,dy)ds 
is a //^”'L'’'^^-martingale. 


While it is important to keep the above setup of martingale problems in mind, we will pass 
on to construct convenient random processes (on a possibly different probability space) for 
ease in calculations. Given that —)• i/ and —)• y, by the Skorohod representa¬ 
tion theorem, we can construct random processes and {Z,X) such that 

£(^K.+i)^^K.+i)) ^ C{Z,X) = 12, and 

(^K.+i)^^K.+i)) ^ (^,^) 

Independently of we construct processes (Z^”L\X^"'L^) and (Z,X) with 

£(zK\xK^) = £(Z,X) = y 

such that X^"'L^) —)• (X, X) a.s. The latter processes are used in this section only to 

shorten certain expressions and write them in more convenient forms to aid calculations. In 
terms of we are able to write the requirement (ii) in the statement of the 

martingale problem as follows: Fix any finite integer r. For any 0 <si<---<Sr<s<t<T, 
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and any choice of bounded Fsi functions 51^ for alH = 1,2, • • • , r, we require 


E 




f f { 4 >{z’"''‘*'^)+a{Z, 

Js Juo 


) - m 


(rifc.+l) ^Jnk.+i) 




) - / V,cj){z, 


{nk,+l) (rife.+l)^ ^{rik.+l) 
u ? -^u 


)-Zu 


du 


c\ 

) T\Zu 1 Zu )} 


(T{\zi ^ -Vu\)l 3 {\Xu ^ 




= 0 


( 6 . 1 ) 


By letting j —>■ 00, we know that 


- ct>{Zt,Xt)\ ^ 0 ( 6 . 2 ) 


if f € H where H := {a G [ 0 ,T] : z^(A(Za,Xa) / ( 0 , 0 )) = 0 }. A similar statement holds when 
t is replaced by s provided that s £ D. The statement that 


E\jy^(^{Z 


(rifc.+i) 

? -^U ) ’ 


(rik.+l) 

)-Zu ^ du- j X,(j){Zu,Xu)- Zudu\^^ ( 6 . 3 ) 


follows from the L^{P) boundedness of supQ<^<7^|Zu ^ \ indexed by j. Hence, our next 

objective is to show that the last term on the left side of ( 16 .ip converges to a limit. 

For (j) G C^{R^ X m 3 ), we will use the notation 


\ 2=1 


3 

W I.—P 1 

2 -^ ' ' 


I WrzzWco 


2-^ dzjdzn 


1/2 






with similar meanings for ||<^j;||oo, ||</'xxlloo, and H^^a-Hoo- 
Fix any function (j) G C'^(M 3 ). On M^^ x E, define 

G{z,x,v,y,^) = {(t){z + a{z,v,^),x) - (l){z,x)}a{\z - u|)/ 3 (|x - y\). 

Claim: 


f\ [ G{zi"^^^"\xi"^^^"\vs,ys,C)f^^"'^2{dv,dy)dQ 
Jo JUo 


lUo 


G{Zs,Xs,Vs,ys,Odidv, dy)dQ\ds 


converges to 0 almost surely. 
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Proof. Denoting the above expression by \ Zg, Xg, 




< D{Zr^'^\xr^ '^',Zg,Xg,fi^^'^^>,^i^^'^^>) + D{Zg,Xg,Zg,Xg,fi^^'^^\fi) 

= h+l2 


for short. Then, /i < Ji + J2 where 




h«fc,+l)x 




io Ji 7 o 

- {(/-(Z, + a{Zg,Vg, 0 ,Xg) - (P{Zg,Xg)}] 

a{\zl ^ ^ - yg\)dQfj,^'^'"^\dv,dii)\ds, and 


J 2 = / I / { 0 (Z, + a(Z„r;„O,^.)--/>(^.,^.)} 
Jo JUo 

dnk,+^) IN ^/, 


{a(|Zr' - Vg\)Pi\Xg^ - 2 /d) - ct{\Zs - Vg\)Pi\Xg - 2 /d)}dQd/iK 

In order to bound Ji, we bound a and /? by 1 in the above expression, and then, we break up 
the rest of the integrand appearing in Ji as follows, and estimate each part separately. 


Ji< f [ [\mzi""^^'Wa{z 

Jo Jn^J r 


,Vg,i),Xg )-(t>{Z: 


^{n^..+l) (nfe.+l) 

, A.q 


')} 


- {cf{Zg + a{zi"'^^^"\vg,0,xi"^^^"^) - cf{Zg,xi"'^^^"^)}] 


T [{f^iZs + oi{Z, 


(nfe +1) 


Jnk.+l) 


(nfe + 1 ) 


Vg, 0 ,X'g ^ ')-cP{Zg,X; ^ ')} 


- {cP{Zg + a{Zg,Vg,0,xi"’‘^^"^) - ^{Zg,xi"’‘^^"^)}] 
+ mZg + a{Zg,Vg, 0 ,xl"'^^^^^) - f){Zg,xl"'^^^^^)} 

- {(fiZg + a{Zg,Vg,(,),Xg) - (j){Zg, Xg)}]\dQdn^'''"o'> d.s 


< 


[ f [\\(t^'zz\\oo\a{z[ \vg,i)\\Z, 

Jo JUn 


(nj. + 1 ) 


I ,„(%■+!) 

\(^[.Zg 


10 JUo 

+ IIV"2||CX)|'-0.^S 

+ \\4>zx\\oo\a{Zs,Vg,f,)\\Xg'°'^ - Xg\]dQn^'^*‘i\dv,dy)ds 
which tends to zero as j —>■ 00, by restricting all processes to Bji as in the previous section. 


Vg,^) - a{Zg,Vg,^)\ 

nk,+i- 
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and using the bounded convergence theorem. Now, we consider J2. 

J 2 = f [ \mZs + a{Zs,Vs, 0 ,Xs)-^{Zs,Xs)} 

Jo Juo 

-Vs\)/ 3 {\xh^'^ -ys\) 

- a{\Zs - Vs\)P{\Xs - ys\)}\dQdy^^'‘^'^\ds 

, Xnki+i) 



< I I U'jUa{Zs,Vs,0 \{\Zs^ -Zs\ + \Xs^ - Xs\}dQdy^^'‘iUs 
lo JUo 

/■* 

<Ci 

lo Jn'^ 


f [ i\z,\ + \vs\){\zl"^^^^'’ - Zs\ + - X,|}d//K)ds 

JO JB 5 


with Cl as a suitable constant, and the last expression above —)■ 0 as j —)• 00 using arguments 
as before. 


Next, we consider I2 where 


h = f\ [ mZs + a{Zs,Vs,0,Xs) - HZs,X,)} 
Jo JUo 

a{\Zs - U5l)/3(|Xs - ys\)dQ{dy^'^'"i^ - dy)\ds. 


We write the expression within absolute value as 

Ejj{<PiZ, + a{Zs,Z^"'^^\i),Xs) - <P{Zs,Xs)}<j{\Zs - 

- {cl>{Zs + a{Zs,Zs,0,Xs) - ^{Zs,Xs)}a{\Zs - Z.|)/3(|X, - X,|)]dQ. 

~ {P'k ■ ) ~ (^fc - ) 

where E refers to expectation with respect to the random variables Zg ^ Xg ^ ,Zs, and Xg . 
Using this, 

h < rE[cii</.;iioo{i^l"'''^ - ^.1 + ii<^iioo(i^.i + \Zg\) 

Jo 

X {\Zg"^^^ - Zg\ + - X,|)}]ds 

which goes to zero as j —)• 00 boundedly by restricting the processes appearing in the above 
expression to a compact set. By bounded convergence theorem, I2 0 . This finishes the proof 
of Claim. □ 
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Next, we observe that 


E sup 
0<t<T 


f j 

Jo JUo 


(nk^+i) 


'Uo 


+ 0({Z. 




) 


(n-ki+i) 


) - m 


(nk^+l) (nfc.+l) 


^)} 


a{\Zs^^ - Vs\)Pi\Xs - ys\)dQd^^"''"j^\‘^ds 


< m\i. 


[ Ed 

Jd2 


z, 


('^fcj+1) 1 2 


+ \vs\'^)dfi^"''°i^ds 


< C 


<Ci 


for a suitable constant C since E [ sup [|zi + \zi 

\o<t<T J 

for a constant Ci which is independent of Uk^ ■ 

E[f\[ G{zi"'‘^^"\xi"'‘^^"\vs,ys,Ofi^"'‘^\dv,dy)dQ 
Jo JUo 


lUo 


G{Zs,Xs,Vs,ys,Od-idv,dy)dQ\ds] 0 


(6.4) 


By (I6.2D - (I6.4h . we conclude that the martingale problem posed by the {uk^ + l)st stochastic 
differential equation converges to a solution of the martingale problem posed by the stochastic 
system 


Zf = z, 


^IL a{Z^,Vs,0\o, ai\zjt-v4)0m-ys\)]ir)dd^^ 

+ [ a{Z^,Vs,0l3{\X^ - ys\)y{dydv)Q{d^)ds 
Juo 

Xf = Xo+ [ Z^ds,. 

Jo 


(6.5) 

( 6 . 6 ) 


We know from Section 3 that the above system, with /r given, has a unique solution. Since 
converges in law to u, and the law of {Z^,X^) is given by u. By the result in Section 
5, we know that if converges in law, then also converges in law to the same 

limit. Hence u = y, and y is a solution of the martingale problem posed by (j6.5|) and (16.61) . 
The proof of Theorem 12.11 is thus completed. 


7 Uniqueness of Solutions 


In this section, we will study the uniqueness of solutions to the martingale problem posed by 
the Enskog equation. Uniqueness is taken in the sense of uniqueness of time marginals of the 
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distribution of the processes (X,Z). We consider the Enskog equation under the additional 
hypothesis that the law of any weak solution of the limit equation admits at each time point 
t, a density with respect to the Lebesgue measure on M®. Such a hypothesis can likely be 
replaced by imposing certain conditions on the functions a and a that are expedient though 
it would take us far away from the physics of the problem. 

Theorem 7.1. Let the hypotheses used in Theorem A2. 1\ hold. In addition, let a be in 
Then, for any fixed t in the interval [0,T], the t-marginal distribution of weak solutions of the 
limit equation (given below) is unique within the class of Borel probability measures on that 
are absolutely continuous with respect to the Lebesgue measure on M®. 


Jo 

+ [ [ a{Zs,Vs,f)a{\Zs-Vs\)l3{\Xs-ys\)i^{dy,dv)Q{d^)ds 

Jo Juo 


Xt = Xo + 



Zgds 


(7.1) 

(7.2) 


wherein the law of {X, Z) is given by n, and the compensator of N is given by ix{dy, dv)Q{dSfidt. 


Proof. Let h be in the space L := x M^). By the Ito formula, we have 


]E[/i(Xi, Zt)] = E[/i(Xo, Zo)] + E / V,/i(X„ Zs) • Zsds 


+ 


¥.[f f {h{Xs,Zs + a{Zs,Vs,0)-h{Xs,Zs)} 

Jo Jd2xH 

a{\Zs - U5l)/3(|X5 - ys\)y.{dy,dv)Q{df,)ds]. 


(7.3) 


If there exists another solution of dZU) and (1721), let us denote it as {X',Z') with its law 
denoted by p. For each s G [0,r], the measures ps and ps denote the s-marginals of p and p 
respectively. 

Using the fact that both {Xt, Zt), {Xfi Zj-) are solutions of (|7.ip and (|7.2I) . and using E to 
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denote expectation of functions of {Xt, Zt) as well as {X^, Z[), we have 


m{Xt,Zt)-m{x[,z[) 


= E 



\y xh{Xsi Zg) ■ Zg 


V^h{X'g,Z',)-Z'g] ds 


+ 


E[ 



{h{Xg,Zg + a{Zg,Vg,i)) - h{Xg,Zg)]a{\Zg - 
ld{\Xg - yg\)n{dy,dv)Q{di)ds\ 

+a(Z',7;;,0) -- ^sl) 

- v's\)p{dy\dv')Q{di)ds\ 


where 


= E 


f [X,h{Xg,Zg) • Z, - V,/i(X', ZO • z'] ds 

Jo 




X,h{X',Z'g)-Z'] ds 


Di := 



10 Jn^ 
I3i\xs 



IQ JO^ JO^ 

IJ{\x'g 


' [h{zg + a{zg,Vg,C))h{zg)]a{\zg - 
yg\)y{dy, dv)y{dx, dz)Q{d^)ds 
[h{Zg + a{Zg,Vg,^)) - h{zg)]a{\zg 
y's\)p{dy', dv')p{dx',dz')Q{di)ds. 


-^^1) 


(7.4) 


(7.5) 


We will first obtain an upper bound for Dq. Define /(x, z) := Xxh{x, z) ■ z. 

Fix any R > 0. If \z\ < R, then / is in the space L. By existence of the second moments of 
supo<t<T \^t\ and supQ<^<ji |Z^'|, we obtain that for any hxed t G [0,r], given any e > 0, we 
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can find a compact set K in M® such that 

\f{x, z)\{iJ.t{dx, dz) + pt{dx, dz)] < ^ 


Therefore, we have 


e/2+ / I / f{x,z){p.t{dx,dz) - pt{dx,dz)}\dt 
Jo Jk 


\Do\ < e 


Let R be large enough so that K C Bfi where Bji denotes the i?-ball in M®. The function 
f{x,z)lK{x,z) can be approximated by a function fs in L by replacing Ik by a smooth 
non-negative function which is identically equal to 1 on iL and whose support is in the 6 
neighborhood of K, 6 > 0, which is contained in B^. As 6 tends to 0, fs —>• fix, pointwise. 
By the bounded convergence theorem, 

lim / / \f{x,z) - fs{x,z)\{pt{dx,dz) + pt{dx,dz)}dt = 0. 

<5^0 70 Jk 

Hence, for all 5 small enough, 

n \f{x, z) - fs{x, z)\{pt{dx, dz) + pt{dx, dz)}dt < e/2. 


The above estimates allow us to conclude that 

ft 

Do <e + 


j I j fs{x,z){pt{dx,dz) - pt{dx,dz)}\dt. 
Since fs £ L, the above inequality ipso facto yields 

|L>o|<e+ / sup I / cl){x, z){pt{dx,dz) — pt{dx,dz)}\dt. 

Jo 4>eL J 


Next, we proceed to estimate Di. On the right side of equation (j7.5p . we add and subtract 

n / / [h{xs,Zs +a{zs,Vs,C)) - h{xs,Zs)\(T{\zs - Vs\)l3{\xs - ys\) 

)2 J’E 

p{dy, dv)p{dx, dz)Q{df)ds 

This enables us to split Di and write it as Gi + G 2 , where 

•= / / / I [h{xs,Zs +a{zs,Vs,C)) - h{xs,Zs)]a{\zs - Vs\)P{\xs - ys\) 

Jo Jn^ Jn^ Je, 

p{dx, dz)p{dy, dv)Q{dffjds 

n / / [h{xs,Zs + a{zs,Vs,i)) - h{xs,Zs)]a{\zs - Vs\)^{\xs - ys\) 

)2 Jid 2 JS 

p{dy, dv)p{dx, dz)Q{dff)ds. 
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and 


G2.= / / [h{xs,Zs + a{zs,v'^,C))-h{xs,Zs)]a{\zs-v'^\)/3{\xs-y's\) 

Jo Jo^ Jo^ J'E 

n{dy', dv')p{dx, dz)Q{d^)ds 

~[ [ [ [iHx's,z's + a{z'„v'„C))-h{x'„z',)]a{\z',-v',\)/d{\x'^-yX 

Jo Jd2 Jd2 Je 

p{dy', dv')p{dx', dz')Q{d^)ds. 

To find an upper bound for Gi, first introduce, for each s G [0,r], the function fg defined on 
by 

fs{x, z)= [h{x, z + a(z, V, 0) - h{x, z)]a{\z - u|)/3(|x - y\)Q{dC)ps{dy, dv). 

JR6xS 


Then Gi can be written as 


n fs{xs,Zs){p{dx,dz) - p{dx,dz)). 
)2 


We will now use the assumption in the statement of the theorem that for each s, ps « A, 
and ps « A, where A is the Lebesgue measure on M®. Hence the product measures ps x Ps 
and Ps X Ps are both absolutely continuous with respect to the Lebesgue measure on 
Therefore, given e > 0, there exists a ^-neighborhood G := Gs of the diagonal in M® X M® such 
that 

{ps X Ps T Ps X ps'){G') < e 
For a suitable constant C > 0, we have 


\Gi\<Ge+ / lG<^{x,z,y,v)[h{x,z + a{z,v,^)) - h{x,z)]a{\z - v\) 

JK6 Jr^xE 


p{\x - y\)Q{d^)ps{dy,dv){ps{dx,dz) - ps{dx,dz)). 


(7.6) 


We approximate the function Iqc by a smooth bounded function (p which takes the value 1 in 
G'^, and zero in G 5 / 2 - With such a choice of (/>, we have 

/ \{^G<^{y,v) - 4>}{x,z,y,v)[h{x,z + a{z,v,i)) - h{x,z)] 

jri2xe: 

sdz - v\)l3{\x - y\)Q{di)\ps{dy,dv){ps{dx,dz) + p{dx,dz)] 

< Ge. (7.7) 

By (17.61) and (17.7p . it follows that 

|Gi|<2Ge-|-| j I (p{x,z,y,v)[h{x,z + a{z,v,^)) - h{x,z)] (7.8) 

7r6 Jr^xE 

a{\z - v\)l3{\x - y\)Q{dP,)\ps{dy,dv){ps{dx,dz) - ps{dx,dz))\. (7.9) 
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The integral 



x,z,y,v)[h{x,z + a{z,v,^)) - h{x,z)]a{\z 


v\)l3{\x - y\)Q{d^)\ns{dy, dv) 


is smooth and bonnded as a fnnction of x, 2 , and therefore it is in L. Hence, (uni) yields 


IGil < 2Ce + 


sup I / 'il:{x,z){yLt{dx,dz) — pt{dx,dz)}\dt. (7-10) 

PeL J 




To bound G 2 , we introduce the function Qs on by 

Qsijj, v) := j[h{x, z + a{z, v, 0) - h{x, z)]a{\z - u|)/3(|x - y\)Q{di)ps{dx, dz). 

Following the ideas used in finding an upper bound for Gi, one can, mutatis mutandis, bound 
G 2 , and conclude that 

|G2|<2Ge+ [ sup| [ 'ijj{y,v){pt{dy,dv) - pt{dy,dv)}\dt. (7.11) 

Jo pGL J 

In summary, we have for any e > 0, 

\E[h{Xt,Zt)-h{XlZ[)\\<hCe + 2, T sup |E[,/.(X„ Z,)] - E[</.(X', Z')]| ds (7.12) 

Jo <t>&L 

for all h £ L, for a suitable constant G > 1. With e being arbitrary, let e —)• 0. By the Gronwall 
Lemma, the uniqueness of the time-marginal distributions of weak solutions is obtained. □ 

Remark 7.1. In the context of martingale problems, uniqueness of distribution of the time 
marginals of solutions implies uniqueness of the law on the path space (JEEI, page 69). Hence, 
Theorem \7.1\ gives us uniqueness in law of the solution of (EH) and (EH on B X D. 


8 Invariant Gaussian density for velocity 


Let {Xs,Zs}sg]R+ be a solution of (l2.7D . (l2.8h . which corresponds to the Enskog equation in 
the kinetic theory of gases. Let MVN (0,/) denote the standard normal distribution on 
where 0 stands for the mean vector, and /, for the 3x3 identity matrix for the variance. In 
the following ’’density of measures” shall be understood relative to the underlying Lebesgue 
measure. 

Theorem 8.1. Let us assume that the law of the initial velocity Zq and that of the initial 
location Xq of EH,EH are independent. Let Zq have MVN (0, 1) distribution. Assume 
that the distribution r]Q of the initial location Xq has density h{x), x G M^. Then the joint 
distribution p{dx,dz) of {Xg, Zs}s£R+ has for all t > 0 density pt{x,y):= ht{x)g{y), where 
g{y) is the density of the normal distribution MVN {0,1), while ht{x) is the density of Xf. 
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Remark 8.1. In particular the marginal velocity Zt at time t is distributed according to the 
MVN (0,/) distribution for all t >0, and is independent of the location Xt for all t >0. 

Proof. Our method of proof relies on guessing the solution and proving that it is indeed the 
solution. We take 


with 


and 


fit{dz, dx) = fj.t{dz 1 

x)rjt{x). 

(8.1) 

p,t{dz x) := g{z)dz 

Vt > 0 

(8.2) 

rjt{dx) := ht{x)dx 

Vt > 0, 

(8.3) 


where ht{x) is a probability density function on which is ascertained below. We will then 
prove that fJ.t{dz, dx) is the distribution of a process ^s}sgr+ which solves (I2.7p . (l2.8p with 
Zg having the same distribution as Zq and solving (12.7p . It would then follow that j^Zgds is 
a Gaussian random variable for all f > 0, and Xt has therefore a density function denoted by 
ht{x). 

Consider := E for any A G and t > 0. It is enough to prove that f/tW = 

00(A) for all A. Using the ltd formula and taking expectation, one obtains 

0,(A) = 0o(A) + r / Li{\^s+^{^s,Vs,i)) _ gi(A,z.)| 

00 iDxDxH ^ 

a{\zs - Us|)/3(|x - y\)pLs{dz, dx)p,s{dv, dy)Q{df)ds. (8.4) 


This is an equation that is satisfied by the characteristic function of Zt where Z is a solution 
of (j2.7p . If iit{dz,dx) is as specified above, then we can write 

0,(A) = 0o(A) + f f _ e*(A,^U| 

Jo 0R6xM6xH ^ 

a{\zs - Vs\)f3{\x - y\)g{z)hs{x)dzdxg{v)hs{y)dvdyQ{df,)ds (8.5) 

which we write as 0o(A) + I- Let us write I as 


I = 


f f 

Jo JrSxS 


s, X, y, f)hs{x)hs{y)Q{dC)dxdyds 


where 


(t>{s,x,y,Ci 


J(\,Zs+a{zs,Vs,^)) 



y\)g{z)dzg{v)dv. 
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Then (f){s, x, y, is 


f ^^i{X,z,+a{zs,Vs,i)) 

JK^xRS 


Vs\)Pi\x - y\)ysidz \ x)ys{dv \ y) 


= f V(N: -t>:i)exp{-l/2(|2,l" + |!..l"))<i^<it. 

{27ry/^ Jr3xr3 

= f T''"-’ - V(|2. - o.l)exp{-l/2(|z:p + 

(zvr)'*/ Jr3xR3 

by Proposition [12 continuing, 

= - ^.1) exp{-l/2(l^,|2 + \vs\^)}dzdv 

{271y^/ Jr3xR3 

= - [ -e^^^’^‘^)a{\zs-Vs\)^{\x-y\)g{z)dzg{v)dv (8.6) 

7r3xR3 

by using conservation of energy in (II. 3P and with z* and v* denoting the post-collision velocities 
corresponding to the pre-collision velocities z and v. It follows that I = —I and hence 1 = 0 
so that (/>f(A) = i?io(A) for all A at all times t > 0. □ 
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